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ABSTRACT 

' Monochromatic opacities from the Opacity Project (OP) (Seaton et al.) have been 

augmented by hitherto missing inner-shell contributions (Badnell & Seaton). OP 
Rosseland-mean opacities, kr, are compared with results from OPAL (Iglesias & 
Rogers) for the elements H, He, C, O, S and Fe. The OPAL data are obtained from 
the website www-phys.llnl.gov/Research/OPAL/index.html. 

Agreement for H is close everywhere except for the region of log(T) ~ 6 and 
| log(i?) — — 1 (R = p/T§ where p is mass-density in g cm" 3 and Tq = 10" 6 x T with T 

y-^- ■ in K). In that region kr(OPAL) is larger than kr(OP) by up to 13%. The differences 

f^S are due to different equations of state (EOS). In the region concerned OP has the H 

ground state undergoing dissolution, leading to a small H-neutral ionization fraction, 
■ while OPAL has larger values for that fraction. A similar difference occurs for He 

at log(i?) ~ —1 and log(T) ~ 6.4, where OP has the He + ground-state undergoing 
dissolution. 

The OPAL website does not provide single-element Rosseland means for elements 
other than H and He. Comparisons between OP and OPAL are made for mixtures 
with X = 0.9, Z = 0.1 and Z containing pure C, O or S. There are some differences: 
at the lower temperatures, say log(T) ^ 5.5, due to differences in atomic data, with 
the OP R-matrix data probably being the more accurate; and at higher temperatures 
mainly due to differences in level populations resulting from the use of different EOS 
theories. 

In the original OP work, R-matrix data for iron were supplemented by data ob- 
tained using the configuration-interaction (CI) code superstructure. The experi- 
ment is made of replacing much of the original iron data with new data from the CI 
code AUTOSTRUCTURE. Inclusion of intercombination lines gives an increase in kr of 
up to 18%. 

The OPAL website does not allow for Z containing pure iron. Comparisons are 
made for an iron-rich mixture, X — 0.9, Z = 0.1 and Z containing C and Fe with 
C:Fe=2:l by number fraction. There are some differences between OP and OPAL for 
that case: the OP 'Z-bump' in is shifted to slightly higher temperatures, compared 
to OPAL. 

Overall, there is good agreement between OP and OPAL Rosseland-mean opacities 
for the 6-elements, but there are some differences. Recent work (Bahcall et al.) has 
shown that helioseismology measurements give a very accurate value for the depth of 
the solar convection zone, i?cz, and that solar models give agreement with that value 
only if opacities at i?cz are about 7% larger than OPAL values. For the 6-element 
mix at Rqz we obtain kr(OP) to be larger than kr(OPAL) by 5%. 

Key words: atomic process - radiative transfer - stars: interiors. 
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1 INTRODUCTION 

Rosseland-mean opacities from the Opacity Project (OP), 
as originally presented in pQ, were in good agreement 



with those from the OPAL project [5] over much of the 
temperature-density plane, but were smaller than those 
from OPAL at high temperatures and densities. Iglesias and 
Rogers 3 offered the explanation that OP was missing some 
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data for inner-shell transitions and that was confirmed in a 
recent paper by the present authors fl]. The results of 
and [I] were for a mixture of 6 elements, H, He, C, O, S and 
Fe with abundances, by number-fractions, given in Table 1: 
we refer to that as the 6-element mix; mass-fractions are 
X = 0.7 for H, Y = 0.28 for He, Z = 0.02 for 'metals'. 

Fig. 1 shows the level of agreement between OP and 
OPAL, for that mix, as obtained in FjQ: values of log(ftR), 
where kr is Rosseland-mean opacity in cm 2 g _1 , are plotted 
against log(T) for 5 different values of log(iZ) where R = 
p/Ti, p = density in g cm" 3 and T 6 = 10 6 x T with T 
in K (fig. 1 of [J shows the behaviour of \og(R) for a few 
typical stellar models). The OPAL data are obtained from 
the OPAL website and the OP data are with inclusion 
of the inner-shell contributions discussed in gj. It is seen 
that the agreement between OP and OPAL is fairly good in 
all cases but that there are some differences. The purpose 
of the present paper is to make more detailed comparisons 
of OP and OPAL for the six elements H, He, C, O, S and 
Fe. Work on the inclusion of OP inner-shell data for other 
elements is in progress. 



2 ROSSELAND MEANS 

Let (Tfc(u) be the cross-section for absorption or scatter- 
ing of radiation by element k, where u = hv /(ksT), v 
is frequency and fee is Boltzmann's constant. For a mix- 
ture of elements with number fractions ft, ^2 k fk = 1, put 
a ( u ) — X/fe/ fe<Tfe ( M ) ( see section f4.ll . The Rosseland-mean 
cross section is <tr where, 
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The Rosseland-mean opacity per unit mass is kr = an./ p, 
where p is mean molecular weight. 



3 EQUATIONS OF STATE 

The populations of the energy levels providing absorption 
or scattering of radiation are determined by the equation of 
state, EOS. OP and OPAL have very different treatments 
of the problem: the former use what is referred to as the 
'chemical picture'; the latter the 'physical picture'. 



3.1 The OP EOS 

Let i be equal to the number of bound electrons in ionization 
stage i. In the OP work the internal partition function for 
stage i is taken to be 6 
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where j specifies an energy level, Qij is a statistical weight, 
Wij an occupation probability and Eij is the total energy of 
level ij. Let 0; be the fraction of atoms in stage i. In the 
OP work (0 and 0) 



where 



. 2-wh 2 / 



3/2 



(4) 



(5) 



and m is the electron mass. Allowance for electron degener- 
acy, and other refinements, is discussed in [Sj and in section 
IV(e) of 0. 

In the OP work the Wij are calculated using methods 
described in 6 which are plausible but not rigorous (some 
modifications of the treatment in [Hj are discussed in [J]). 
The OP Wij become small for levels which are sufficiently 
highly excited (referred to as 'level dissolution'), which en- 
sures convergence of the summation in (T3J. In many cir- 
cumstances the Boltzmann factors in (T3J, exp[— Eij/(k&T)\, 
become small long before there is a cut-off due to the Wij 
becoming small: in those circumstances the calculated opac- 



3.2 The OPAL EOS 

The OPAL approach to the EOS problems is based on the 
many-body quantum statistical mechanics of partially ion- 
ized plasmas (see |33| and |35| for a summary of later work) . 
The level populations obtained by OPAL may be expressed 
in terms of the W factors of OP. Comparisons of results from 
OPAL and OP are given in I21| for H and H + and in [2] for 
hydrogenic C. It is found that, compared with OP, OPAL 
has larger populations in the more highly excited states. 

The case of hydrogenic C was discussed further in pP 
where it was found that OPAL gave W factors to be surpris- 
ingly large for states which had mean volumes larger than 
the mean volumes available per particle in the plasma. 



4 OPACITIES FOR MIXTURES 

In principle, the level populations for any one chemical ele- 
ment depend on the abundances of all other elements present 
in a plasma. However, it would not be practicable to make 
complete ab initio calculations of opacities for every mix- 
ture which may be of interest. Some approximations must 
therefore be made. 



4.1 The OP approach 

The OP Wij depend on the ion micro-field which should, 
of course, depend on the chemical mixture. The approxima- 
tion is made of using a micro-field independent of mixture 
(in practice, that for fully-ionized H and He with X — 0.7, 
Y — 0.3). Monochromatic opacities are then calculated for 
each chemical element, as functions of frequency, on a grid 
of values of T and N e (in practice, usually with intervals of 
<Hog(T) = 0.05 and 51og(iV c ) = 0.5). The monochromatic 
opacities are then simply added together for the calcula- 
tion of Rosseland means for mixtures. In most cases that 
procedure does not lead to significant error, but there are 
exceptions (see, for example, Section I7.3I - ). 
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4.2 The OPAL approach 

The OPAL work makes use of interpolation procedures 
based on the concept of 'corresponding states' (see 
Rosseland-mean opacities are available from the OPAL web- 
site for 'metal' mass-fractions Z ^ 0.1. In general the use 
of interpolations does not introduce any important errors 
(see but there may be some significant errors near the 
edges of the domains of the tables provided. 



5 ATOMIC PHYSICS 
5.1 R-matrix calculations 

Most of the atomic data used in the original OP work, were 
obtained using the R-matrix method, RM |10|. For a system 
containing N electrons we use the wave-function expansion 

* = A ^ ipp (xi , ■ ■ ■ ,xjy_i) x 9 p (x N ) 
v 

xi, ■ ■ ■ ,Xiv) x c m (6) 

m 

where: A is an anti-symmetrisation operator; Xi is a space- 
and-spin co-ordinate for electron i; ip p is a function for a 
state k of the (N — l)-electron core (usually calculated us- 
ing a configuration-interaction code); 8 P is an orbital func- 
tion for an added electron; $ m is a function of bound-state 
type for the iV-electron system; c m is a coefficient. In the 
R-matrix method the functions 8k and coefficients c m are 
fully optimised. 

Further details on the R-matrix calculations are given 
in [TT] and E]. 



5.2 SUPERSTRUCTURE and AUTOSTRUCTURE 

The configuration interaction (CI) codes SUPERSTRUCTURE 
(SS) 23] and autostructure (AS) use expansions 

* = £>™ (7) 

m 

where the $ m are one-configuration functions. The code 
SS is used only for calculating energies of bound states 
and radiative transition probabilities: AS, which was devel- 
oped from SS, includes calculations for auto-ionizing states, 
auto-ionization probabilities, and cross-sections for photo- 
ionization. 

In using the CI codes for a level of a given configuration 
C, we usually attempt to include in the expansion (6) at 
least all of the states belonging to the complex of which C 
is a member. 



5.3 Methods used in OPAL 

OPAL uses single-configuration wave functions, 

* = (8) 

with one-electron orbitals calculated using potentials ad- 
justed empirically such as to give best agreement with ex- 
perimental energy-level data |16l I17| . 



6 SPECTRUM LINES 

Inclusion of contributions from large numbers of spectrum 
lines is of crucial importance. 

6.1 Inclusion of lines 

Our procedure is, first, to calculate the monochromatic opac- 
ity cross-section a(y) including only continuum processes 
(photo- ionization and -detachment, free- free processes and 
scattering) and to determine a mean back-ground <tb(i / ) 
such that, at each frequency point v, a(y) ^ <tb(^). Let 
<Ji,{y) be the contribution from a spectral line. We define a 
quantity testl and, at each frequency mesh-point, include 
the line contribution if 

> testl x ctb (y) - (9) 

We take testl to be small (usually testl= 10~ 8 ): the con- 
tribution from each line may be small but the number of 
lines may be very large. 

6.2 Fine structure 

Inclusion of the fine structures of the spectrum lines can lead 
to a re-distribution of oscillator strength which can give a 
significant increase in or defined by Q . In the original OP 
work fine-structure was allowed for by methods described 
in section 4.5 of 0: simple LSJ-coupling formulae were used 
together with empirical estimates of averaged spin-orbit pa- 
rameters so as to determine level splittings. It was checked 
that the method was very stable, in that variation of the 
parameters over a wide range gave hardly any change in the 
calculated values of or. However, the method did not allow 
for the inclusion of intercombination (spin-forbidden) lines 
which can be important for highly ionized systems. 

6.3 Line profiles 

The quantity or can be sensitive to the widths of the spec- 
trum lines, which are due to: radiation damping; thermal 
Doppler effects; and pressure broadening. Both OP and 
OPAL use similar empirical formulae for the pressure widths 
based on a theory originally proposed by Griem 1181 . They 
differ in that OPAL adopts parameters determined from 
available experimental data |19| while OP uses results from 
ab-initio calculations |20l 1211 122| . In the fairly small num- 
ber of cases for which direct comparisons have been made 
|21| the agreement between OP and OPAL pressure-widths 
is very close. 

It may be noted that or is sensitive to pressure broad- 
ening only for intermediate values of the density: at suffi- 
ciently low densities the pressure- widths are small compared 
with those due to radiation damping and Doppler effects; 
while at high densities the lines are blended to form a quasi- 
continuum. 



7 HYDROGEN AND HELIUM 

The basic physical data for H and He (energy levels, radia- 
tive transition probabilities, and cross-sections for photo- 
detachment, photo-ionization and free-free transitions) are 
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known accurately and theories of pressure-broadening of 
lines for those elements are well developed and should be 
quite reliable. Fig. 2 gives OP and OPAL results for a H/He 
mixture with mass-factions of X — 0.7 for H and Y — 0.3 
for He. As is to be expected, the agreement between OP and 
OPAL is generally very close but there is a region with some 
differences for log(i?) = — 1 within the vicinity of log(T) = 6. 
We discuss that region further. 

7.1 Hydrogen 

The top part of Fig. 3 shows log(/tR) for pure H with 
log(i?) = — 1 and log(T) = 5 to 7 and the lower part shows 
the first derivative, 

aiog(*H)/0]og(T), (10) 

calculated in the approximation of using first differences. 
The largest difference between OPand OPAL occurs for 
log(T) ~ 6.0 and in that region there is a marked change in 
the behaviour of 91og(/tR,)/c?log(T), which becomes practi- 
cally constant for log(T) > 6.5. There are three main contri- 
butions to hydrogen opacities in the region considered: (1) 
scattering of photons by free electrons; (2) electron-proton 
free-free transitions; (3) hydrogen bound-free transitions. 
Contributions from (3) become very small for log(T) signif- 
icantly larger than 6.0, as H becomes fully ionized. 

OP calculations are first made on a mesh of values of 
(log(T), log(jVe)) and interpolations to required values of 
(log (T), log (p)) are then made using the code OPFIT.F |23|. 
Table 2 gives results for the Hydrogen ionization equilibria 
at log(T) = 6.0: values of log(ii); values of log(iV c ); 0(1) 
(fraction of neutral H); and W(l) (H ground-state occupa- 
tion probability) from OP. The table also includes values 
of log(iVe) and from OPAL (data kindly supplied by 
Dr C. A. Iglesias). The trends for OP may be described as 
follows: with increasing density, (j>(l) initially increases due 
to the factor N e /U e in (3) (pressure recombination) but for 
the larger values of density W(l) becomes small leading to 
smaller values of <f>(l) (pressure ionization). The values of 
4>(1) from OPAL are much larger than those from OPAL 
and the additional amount of neutral H explains the differ- 
ence between the OP and OPAL opacities. 

It will be noted that at the lower densities of Table 2 the 
values of 4>(1) from OPAL continue to be larger than those 
from OP. The reason is that the OPAL values of the W(n) 
for n > 1 are larger than the OP values. However, those 
differences in ionization equilibria at the lower densities do 
not have much effect on the Rosseland-mean opacities. 

7.2 Helium 

Fig. 4, for Helium, shows results similar to those of Fig. 3, for 
Hydrogen. At log(i?) = —1 the maximum difference between 
the OP and OPAL opacities occurs for log(T) ~ 6.4. That 
is the region where the OP occupation probability for the 
He + ground-state becomes small. The OPAL populations 
for He + will be much larger than those from OP. 

7.3 Smoothness 

A further feature of Figs 3 and 4 may be noted. Both OP and 
OPAL give rise to somewhat irregular appearances for the 



derivatives 91og(KR,)/91og(T). Dr Iglesias informs us that 
when producing a pure H or He table for OPAL it is nec- 
essary to do some interpolating, which might explain the 
irregularities in the OPAL data. The OP work also involves 
interpolations. A check run was made for H using intervals 
of <51og(T) = 0.025 and 51og(jV e ) = 0.25, that is to say one 
half of those normally used, and it was found to give close 
agreement with the results presented in Fig. 3. 

7.4 Summary on H and He 

The differences between OP and OPAL for level populations 
(and, more generally, for equations of state, EOS) have been 
discussed in a number of previous papers (|24|. |1]) and 
will not be discussed further here. Our present concern is 
only with effects on opacities. For pure H and pure He the 
largest differences in opacities, at log(i?) = —1, are 13% for 
H at log(T) = 6.0 and 20% for He at log(T) = 6.4. For 
smaller values of log(i?) the differences are much smaller 
(see Fig. 2). When other elements are included (the 'metals') 
the differences between OP and OPAL are generally much 
smaller, because in the regions concerned the contributions 
to opacities from those other elements will generally be much 
larger than those from H and He. Thus Fig. 1 shows no sign 
of the discrepancy between OP and OPAL shown on Fig. 2. 



8 CARBON AND OXYGEN 

For the C and O RM calculations we include in ||SJ all states 
ip p which belong to the ground-complex of the (N — 1)- 
electron system, that is to say all states with the same set of 
principal quantum numbers as the ground state. Thus the 
system with N = 6 electrons has a ground configuration of 
ls 2 2s 2 2p 2 and the ground-complex of the (N — l)-electron 
system contains the configurations ls 2 2s 2 2p, ls 2 2s2p 2 and 
ls 2 2p 3 . With those choices for the ip p , all bound states of 
C and O can be represented by expansions of the type JSJ. 
R-matrix calculations were made for all bound states with 
an outer electron having quantum numbers n£ with n ^ 10 
and £ ^ lmax with LMAX varying between 2 and 4: states 
with £ >lmax were treated in hydrogenic approximations. 

Independent evaluations of OP oscillator strengths for 
C, N and O have been made in 1251 . The OP values are 
found to compare favourably with experimental data and 
with data from other refined calculations. For most transi- 
tions in C, N and O the OP data are recommended in |25| as 
the best available. For C and O, in addition to the R-matrix 
data, we now include data for inner-shell transitions from 

m 

We recall that the OPAL website [Sj provides 
Rosseland-mean opacities for 'metal' mass-fractions 
Z ^ 0.1. Figs 5 and 6 compare OP and OPAL opacities 
for H/C and H/O mixtures with X = 0.9 and Z = 0.1. 
Agreement is fairly close. 

For the lower temperatures, log(T) < 5.5, the differ- 
ences between OP and OPAL may be due to differences in 
atomic data, with the OP RM data probably being the more 
accurate. The maxima at higher temperatures (log(T) ~ 6.0 
for C and 6.2 for O) are due to transitions with K-shell ini- 
tial states (electrons with principal quantum number n = 1). 
For the higher densities in those regions, log(-R) > —2, the 
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opacities from OP are larger than those from OPAL, which 
is an EOS effect. 

In the regions concerned the dominant ionisation stages 
for C and O are H-like and He-like and for OP dissolution 
occurs for outer electrons having n > 3, while for OPAL 
it occurs at considerably larger values of n; hence OPAL 
will have more electrons in highly excited states, and less in 
states with n — 1, giving smaller opacities. For the reasons 
given in section 3.2 we consider that OP may give the more 
accurate results for the positions at which dissolution occurs. 



9 SULPHUR 

For ions with numbers of electrons TV ^ 11 the ground- 
complexes for the (TV — l)-electron cores contain only Is, 2s 
and 2p electrons but for TV > 11 3d electrons must also be 
included and the number of configurations in the complexes 
can become very large. For such values of TV that number 
was, in general, too large for the inclusion of all such states in 
the RM expansions and some further approximations were 
required. Thus even for the case of TV = 13 the states 3s 2 , 
3s3p ands 3p 2 were included but 3p3d and 3d 2 were omitted. 
For the case of S, with TV ^ 16, such omissions occurred only 
for the earlier ionization stages where the omitted core states 
would be rather high and their omission would be unlikely 
to lead to serious error. 

Fig. 7 compares OP and OPAL opacities for an H/S mix 
with X = 0.9 and Z — 0.1. The level of accord is similar 
to that observed for C and O. For S, kr(OP) is larger than 
kr(OPAL) at the higher densities in the K-shell region at 
log(T) ~ 6.7. That is an EOS effect. 



10.1 Summary of iron data used 

Calculations have been made in both LS coupling and in 
intermediate coupling (IC). Table 3 gives a summary of the 
atomic data for iron used in various stages of the OP work. 



10.1.1 Data from Kurucz 

For the first few ionization stages, TV = 21 to 26, data from 
Kurucz were used |27| . for nearly 7 million lines. The data 
were computed using the code of Cowan \'2iS\. The contribu- 
tions from those stages is not normally of major importance 
for calculation of the Rosseland means but can be of crucial 
importance for radiative accelerations |29|. 



10.1.2 R-matrix calculations, RM 

RM calculations for energy- levels, oscillator strengths and 
photo-ionization cross-sections were made for all ionization 
stages of iron but for TV = 21 to 26 only the photo-ionization 
data were used (much larger amounts of data being available 
from Kurucz). In Table 3 the numbers of levels and lines are 
given separately for TV = 21 to 28, TV = 14 to 20 (the region 
of the Z-bump) and for TV = 2 to 13. For TV = 2 to 13 the 
data from RM are probably the best available. 



10.1.3 Use of SUPERSTRUCTURE 

The RM method was not capable of providing the very large 
amounts of data in the region of the Z-bump, and further 
data for over 3 million lines were calculated using SS I3()|. 
and used in 0. 



10 IRON 

Iron plays a special role in calculations of opacities for typ- 
ical cosmic mixtures, in consequence of its comparatively 
high abundance and the very large numbers of lines in its 
spectra. As the nuclear charge Z increases all levels in a com- 
plex eventually become degenerate (in the limit of Z —* oo 
the one-electron energies, when relativistic effects are ne- 
glected, depend only on the principal quantum numbers 
n). With increasing Z the levels in a complex move down- 
ward and the number of bound levels in a complex there- 
fore increases. Consider iron with TV =16 (Fe XI) which 
has a ground configuration 3s 2 3p 4 and a ground-complex 
containing all configurations of the type ?>s x ip v 'id z with 
x + y + z — 6. Runs with AS give the number of spectroscopic 
terms which belong to that complex, and which have ener- 
gies below the ionization limit, to be equal to 721. There are 
large numbers of radiative transitions between such states 
(the '3 — > 3' transitions) and from such states (the '3 — > n' 
transitions with n > 3). It was first shown in the OPAL 
work that such transitions in iron ions with TV = 14 to 19 
give rise to an important feature in the Rosseland-mean, at 
log(T) ~ 5.2,which has come to be known as the 'Z-bump' 
(see Fig. 1). 



10.1.4 Inner-shell data from AUTOSTRUCTURE 

More recently, inner-shell data from AS, for iron and other 
elements, have been included In |1] it was found that, 
for inner-shell data, use of IC, in place of LS coupling, did 
not give any significant change in Rosseland means. Here we 
consider only LS coupling inner-shell data. 



10.1.5 Further data from AUTOSTRUCTURE 

In the present paper we consider the replacement of data 
from RM and SS, in the regions of the Z-bump (TV = 14 to 
19), with further outer-shell data from AS. Calculations were 
made using both LS coupling and IC with over 30 million 
lines for the latter case. 



10.1.6 Photo-ionization 

The numbers given under the heading 'PI' in Table 3 cor- 
respond to the numbers of initial states for which photo- 
ionization cross-sections were calculated: each cross-section 
contains data for all available final states. 
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10.2 Use of new AS data in LS coupling 

For ionization stages giving the Z-bump, N — 14 to 19, we 
tried replacing all old (RM plus SS) data with the new AS 
data in LS coupling. We obtained Rosseland means close 
to those from the original OP work which provided a good 
independent check. The only significant differences were in 
fairly restricted regions where photo-ionization had been es- 
timated using SS data (the SS estimates were rather crude 
- see EH). 



10.3 Intercombination lines 

We replaced all old (RM plus SS) data for N = 14 to 19 with 
new AS data in IC. The differences with Rosseland-means 
from the previous subsection (AS data in LS coupling) was 
almost entirely due to inclusion of the intercombination lines 
in AS. Fig. 8 shows the percentage increase in the mean 
which results from inclusion of the intercombination lines 
(up to 18% for log(-R) = —4). The inclusion of those lines 
gives improved agreement with OPAL near the maximum 
of the Z-bump. But, on a log-log scale, our final results for 
the 6-element mix do not appear greatly different to the eye 
from those shown in Fig. 1. 



10.4 Many-electron jumps 

Selection rules state that radiative transitions can occur 
only between configurations which differ in the states of one 
electron. In configuration-interaction calculations the states 
may be given configuration labels but transitions can occur 
between states having configuration labels differing by more 
than one electron (the 'many-electron jumps'). The effect 
can be a further re-distribution of oscillator strength. We 
have made runs in which the many electron jumps are omit- 
ted. Fig. 9 shows the percentage increase which results from 
inclusion of those jumps. The effect is not large (no more 
than 8%). It is not included by OPAL. 



10.5 Results for an iron-rich mixture 

It has already been noted that the OPAL website can- 
not provide opacities for mixtures with Z > 0.1. A further 
restriction is that it cannot provide data for cases in which 
Z is pure iron (that restriction arises because such cases 
can be very sensitive to behaviours near minima in the iron 
monochromatic opacities). About the most iron-rich case for 
which we can obtain OPAL data is X — 0.9, Z — 0.1 and Z 
containing C and Fe in the ratio of 2:1 by number fraction. 

Fig. 10 compares OP with OPAL for the iron-rich mix- 
ture. 

At low densities inclusion of iron intercombination tran- 
sitions in the OP work lead to an increase in kr in the region 
of the maximum of the Z-bump, giving close agreement be- 
tween OP and OPAL in that region. But the position of the 
Z-bump given by OP is shifted to slightly higher tempera- 
tures, compared with OPAL. 

At higher densities there are some differences between 
the OP and OPAL results shown on Fig. 10. At those densi- 
ties there are some suggestions of rather irregular variations 



in the OPAL data which might be due to interpolation er- 
rors near the edges of the domain of validity for the OPAL 
interpolations. 



11 THE SOLAR RADIATIVE INTERIOR 

Computed solar models are sensitive to radiative opacities 
in the solar radiative interior, the region below the base of 
the solar convective zone at Rcz- It is shown in two recent 
papers f|31|. |32|1 that helioseismology provides remarkably 
accurate measures of Rcz, and hence stringent tests of the 
accuracy of solar models. In those papers it is noted that 
recent work leads to revisions in solar element abundances 
and that, when those revisions are taken into account, there 
are significant differences between values of Rcz obtained 
from helioseismology measurements and from solar models 
calculated using OPAL opacities. In |31|. from a study using 
envelope models, it is argued that in order to obtain a correct 
density profile it is necessary to adopt opacities larger than 
those from OPAL by an amount of 19%; while in 1321 . from 
a study of full evolution models, it is found that an increase 
by 7% is required. 

At Rcz the best estimates |32| of the temperature and 
density are log(T) = 6.338 and log(p) = —0.735 giving 
log(ij) = —1.75. In Fig. 11 we show the percentage differ- 
ences, (OP - OPAL), in kr against log(T) for the 6-element 
mix at log(ii) = -1.75. At log(T) = 6.338 we find acr(OP) 
to be larger than Kr(OPAL) by 5.0% which is fairly close 
to the value of 7% required in |32|. We find it to be highly 
unlikely that the OPAL opacities could be in error by as 
much as 19% as suggested in |31|. 



12 SUMMARY 

There are two main steps in opacity calculations. 

(i) The EOS problem, determination of the popula- 
tions of all species in a plasma which can lead to absorption 
of radiation. 

(ii) The atomic physics problem, obtaining the 
atomic data which control the efficiencies of the the radiative 
processes. 

12.1 The EOS problem 

The OPAL and OP approaches to the EOS problem are very 
different (see Section 3). 

We encounter two cases for which differences between 
OP and OPAL opacities are due to different treatments of 
the EOS problem. 

The first is for log(_R) ~ -1 and log(T) ~ 6.0 for H 
and 6.4 for He (see Section 3). For those cases OP gives the 
ground states of H° and He + to be undergoing dissolution, 
leading to pressure ionisation, while OPAL has much larger 
populations for those states. The biggest difference is for He 
at log(T) ~ 6.4 and log(-R) = -1, for which kr(OPAL) is 
larger than kr(OP) by 20%. We are not aware of any ex- 
perimental or observational evidence in favour of one result 
or the other. 

The second case concerns K-shell transitions in 'metals' 
at the higher densities, where OPAL has larger populations 
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in excited states, and hence less in states with n = 1 giving 
smaller opacities. That causes the 5% difference between 
acr(OP) and kr(OPAL) at conditions corresponding to those 
at the base of the solar convection zone (see Section 11). The 
helioseismology evidence tends to favour the results from 
OP. 



12.2 Atomic data 

It is essential to consider both the quantity and the quality 
of the atomic data used. 

The OP work was started, a bit more than 20 years ago, 
with the ambitious intention of computing all of the required 
atomic data using the sophisticated R-matrix method, RM. 
The OPAL work was started, at about the same time, with 
the less ambitious approach of using single-configuration 
wave functions computed with the aid of parametric poten- 
tials adjusted empirically so as to give best agreement with 
experimental energies. The OPAL work quickly led to the 
discovery of the 'Z-bump', a feature at log(T) ~ 5.2 pro- 
duced by a very large number of iron lines. OP was unable 
to compute data for such large numbers of lines using the 
RM method, but was able to make supplementary computa- 
tions using the code superstructure which includes some 
allowance for configuration-interaction effects. The first pub- 
lished OP opacities had substantial differences from OPAL 
only at rather high temperature and densities and in [3] it 
was suggested that OP was missing some data for inner-shell 
transitions. That suggestion is confirmed as correct in our 
recent paper The RM and SS codes are both unsuitable 
for handling the inner shell processes and we therefore used 
the code autostructure, AS. In the present paper we also 
include large amount of outer-shell data from AS. 



13 CONCLUSIONS 

Our main conclusion is that there is good general agreement 
between OP and OPAL opacities for the 6 elements H, He, 
C, O, S and Fe. We expect to obtain equally good agreement 
for other elements once work is complete on the calculation 
of further inner-shell data. 

There are some indications that the OPAL interpolation 
procedures described in |5] may introduce some errors near 
the edges of the domains of validity which are claimed. 

At lower temperatures, log(T) < 5.5, there are some 
modest differences between OP and OPAL which may be a 
consequence of the greater sophistication of the OP atomic- 
data work. At higher temperatures and higher densities, 
log(7?) > —2, the OP opacities tend to be larger than those 
from OPAL in consequence of the use of different equations 
of state. There is some helioseismology evidence that the OP 
results are the more accurate. 
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Table 1. The 6-element mixture of references and 



Element 


Number fraction 


H 


9.071(-1)* 


He 


9.137(-2) 


C 


4.859(-4) 





9.503(-4) 


s 


9.526(-5) 


Fc 


3.632(-5) 



*9.071(-1) = 9.071 x lCT 1 . 



Table 2. Hydrogen ionization equilibrium at log(T) = 6.0. 



OP OPAL 



\og(R) 


log(iVe) 






W(l) 


log(iVe) 






-2.776 


21.00 


5.304(- 


-4) 


0.960 


21.00 


1.090(- 


-2) 


-2.276 


21.50 


1.324(- 


-3) 


0.863 


21.49 


1.506(- 


-2) 


-1.775 


22.00 


2.689(- 


-3) 


0.575 


22.00 


2.101(- 


-2) 


-1.276 


22.50 


2.751(- 


-3) 


0.190 


22.49 


2.546(- 


-2) 


-0.776 


23.00 


1.491(- 


-3) 


0.034 


22.98 


4.164(- 


-2) 


-0.276 


23.50 


6.424(- 


-4) 


0.005 


23.45 


1.069(- 


1) 



Table 3. Atomic data for iron. 





N 


Method 


Coupling 


Levels 


Lines 


PI 


21- 


-26 


Kurucz 


IC 


65 427 


6 920 198 




20 


26 


RM 


LS 


8 925 


368 445 


4 479 


14 


20 


RM 


LS 


6 502 


276 450 


4 639 


2- 


13 


RM 


LS 


4 232 


151 974 


3 559 


14 


19 


SS 


LS 


134 635 


3 295 773 


2 537 


14 


19 


AS, outer 


LS 


371 483 


4 594 253 


8 006 


14 


19 


AS, outer 


IC 


1 011 266 


30 555 846 


20 851 


2- 


14 


AS, inner 


LS 


(7 688)* 


2 145 442 


1 339 



* Lower (initial) levels only. 
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Figure 1. Comparisons of log(KR,) from OP and OPAL for 
the 6-element mixture of Table 1. OP from 4 , OPAL from 
EJ. Curves are labelled by values of log(-R) where R = p/T$, 
p is mass density in g cm" " 3 and T 6 is 10 6 x T with T in K. 




Figure 2. Comparisons of log(fi;R,) from OP and OPAL for 
a H/He mixture with mass fractions X — 0.7 for H and 
Y — 0.3 for He. Curves are labelled by values of log(7?). 
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log(T) 

Figure 3. Upper figure: log(Krt) for pure H, log(_R) = 
— 1, log(T) = 5 to 7. Lower figure: the derivative 
<91og(KR,)/(91og(T) calculated using first differences. 
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Figure 4. Upper figure: log(/vR.) for pure He, log(i?) = 
— 1, log(T) = 5.7 to 7.5. Lower figure: the derivative 
<91og(KR,)/(91og(T) calculated using first differences. 
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Figure 6. Comparisons of log(/tR.) from OP and OPAL for a 
H/O mixture with mass fractions X — 0.9 for H and Z — 0.1 
for O. Curves are labelled by values of log(_R). 
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Figure 8. Percentage increase in kr, for the 6-element mix, 
which results from the inclusion of iron intercombination 
lines. Curves are labelled by values of log(_R). 
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Figure 9. Percentage increase in k r , for the 6-element mix, 
which results from the inclusion of iron many-electron-jumps 
(MEJ). Curves are labelled by values of log(_R). 



X=0.9, Z=0.1. H, CandFe. C:Fe=2:1 




log(T) 

Figure 10. Comparisons of log(/tFt) from OP and OPAL 
for an iron-rich mixture: X = 0.9, Z = 0.1 and C:Fe=2:l by 
number fraction. Curves are labelled by values of log(i?). 
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